Caloric curve of star clusters 



(N 

o 

(N 









Lapo Casetti^'^'LJ and Cesare Nardini^'^-^^LlJ 

^ Dipartimento di Fisica e Astronorma and Centra per lo Studio delle Dinamiche Complesse (CSDC), 

Universttd di Firenze, via G. Sansone 1, 1-50019 Sesto Fiorentino (FI), Italy 

^Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Firenze, 

via G. Sansone 1, 1-50019 Sesto Fiorentino (FI), Italy 

' Laboratoire de Physique, Ecole Normale Supeneure de Lyon, 46, allee d'ltalte, F-69007 Lyon, France 

(Dated: February 25, 2013) 

Self-gravitating systems, like globular clusters or elliptical galaxies, are the prototypes of many- 
body systems with long-range interactions, and should be the natural arena where to test theoretical 
predictions on the statistical behaviour of long-range-interacting systems. Systems of classical self- 
gravitating particles can be studied with the standard tools of equilibrium statistical mechanics, 
provided the potential is regularized at small length scales and the system is confined in a box. The 
confinement condition looks rather unphysical in general, so that it is natural to ask whether what 
we learn with these studies is relevant to real self-gravitating systems. In order to provide a first 
answer to this question we consider a basic, simple, yet effective model of globular clusters, the King 
model. This model describes a self-consistently confined system, without the need of any external 
box, but the stationary state is a non-thermal one. In particular, we consider the King model with 
a short-distance cutoff on the interactions and we discuss how such a cutoff affects the caloric curve, 
i.e. the relation between temperature and energy. We find that the cutoff stabilizes a low-energy 
phase which is absent in the King model without cutoff; the caloric curve of the model with cutoff 
turns out to be very similar to that of previously studied confined and regularized models, but for 
the absence of a high-energy gas-like phase. We briefiy discuss the possible phenomenological as 
well as theoretical implications of these results. 
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I. INTRODUCTION 

Many objects in the universe, such as globular clusters, elliptical galaxies, or even clusters of galaxies, can be 
modeled as systems of classical particles mutually interacting via gravitational forces, as long as other interactions are 
negligible compared to the gravitational ones [ij . From a theoretical point of view, these systems also appear as the 
prototype of systems with long-range interactions, whose statistical properties have received a considerable attention 
in the last years [3, |3[. Indeed, the unscreened nature of the gravitational potential makes it truly long-range in any 
situation and should make self-gravitating systems the ideal candidates where to test theoretical predictions on the 
equilibrium and non-equilibrium statistical behaviour of many-particle systems with long-range interactions, like the 
presence of negative microcanonical specific heats, non-homogeneous equilibrium states, long-lived quasi-stationary 
states [1|. 

However, as far as equilibrium properties are concerned, self-gravitating systems suffer from a serious drawback: 
they do not have "true" thermal equilibrium states, neither in the canonical nor in the microcanonical ensemble. 
This lack of thermal equilibrium is due to two main reasons (see e.g. [l|, lj-0] for more detailed discussions), (i) The 
gravitational potential is singular for vanishing distance between two particles, so that entropy and free energy are 
infinite for a system of more than two particles [8|: physically this means that, at equilibrium, (at least part of) the 
system would collapse in states with infinite density, (ii) If the particles have a thermal (i.e., Maxwellian) velocity 
distribution they will tend to escape the system since the escape velocity is finite unless the mass of the system itself 
is infinite. Equivalently, a thermal velocity distribution in a self-gravitating system in unbounded space is possible 
only if the system has infinite mass ^J . 

From a physical point of view the first problem admits a simple solution. No real system exists where the only 
non-negligible interaction is classical gravity at all length scales: either the interacting "particles" are macroscopic 
bodies like planets or stars or galaxies, or quantum effects must be taken into account below a certain length scale. 
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In both cases, a length scale exists (the size of the bodies or the scale where quantum effects set in) below which the 
potential has no longer the classical gravitational form. If one is not interested in small-scale details, the potential 
can be regularized e.g. by replacing it with a softened one at short distances (see e.g. [3, 0|) or by directly considering 
self-gravitating fermions [W\ . 

The second problem is tougher, and one is (somehow artificially) forced to confine the system in a box so that 
statistical equilibrium exists. On physical grounds such an approximation may be reasonable in those cases where the 
evaporation rate is slow compared to the other time scales in the system [1C|, and this happens in many astrophysical 
systems [l|, |ll|; however, in some situations it may also induce nonphysical features, as we shall see in the following. 
We note that also when computing the entropy or free energy of a fluid or a gas one encloses it in a box. However, 
not only in that case this is precisely what is done in a laboratory — at variance with the case of a galaxy or of a star 
cluster — but, being the interactions short-ranged and the equilibrium states homogeneous, the size, the shape and the 
position of the container have no effect on the bulk properties of the fluid. This is not the case for self-gravitating 
systems, the forces being long-ranged and the equilibrium states non-homogeneous. 

Once regularized and confined, self-gravitating systems do have thermal equilibria, both in the canonical and in the 
microcanonical ensemble [121 |l3| , so that they can be approached with the standard tools of equilibrium statistical 
mechanics, like mean-field theory, study of simplified models, and numerical simulations. In the last decades a lot of 
models, ranging from minimalistic toy models to more realistic ones, have been studied from this point of view (we 
shall discuss some of these results, referring to the original papers, in Sec. |ll| for a review see e.g. @, [Tj]). These 
studies have remarkably shown that the collective behaviour of all these models in the microcanonical ensemble is 
qualitatively almost the same. At high energy there is a gas-like phase; lowering the energy a negative specific heat 
phase, where temperature grows while decreasing the energy, is found; this is the phase dominated by gravity. At even 
lower energies the specific heat becomes positive again and remains positive until the minimum energy is eventually 
reached. In this low-energy phase the physics is dominated by the regularization. Also the canonical behaviour of 
these systems is essentially the same and is much simpler than the microcanonical one: there is a discontinuous phase 
transition from the high-temperature gas-like phase to the collapsed phase dominated by regularization. An interval 
in energy, containing the region where the specific heat is negative, is avoided in the canonical ensemble. This is a 
clear manifestation of ensemble inequivalence. 

It is important to stress that even if equilibrium states exist, nothing guarantees that they are actually reached by 
the dynamical evolution of self-gravitating systems regardless of the initial conditions: indeed, long-range interacting 
systems have long-lived off-equilibrium states (called quasi-stationary states, QSSs) where the dynamics can remain 
trapped for times rapidly growing with the number of degrees of freedom [^ so that relaxation to thermal equilibrium 
(if any) may take extremely long times (for a study of this phenomenon in a toy model of a gravitating system 
see [l5|). This, in addition to the non-physical assumption of confinement, might suggest the study of equilibrium 
properties of self-gravitating systems has only a purely academic interest as a mathematical problem. However, real 
self-gravitating systems such as globular clusters or elliptical galaxies do show remarkably similar collective properties 
[31 suggesting that some kind of relaxation has occurred. This makes at least reasonable, although not firmly and 
rigorously grounded, a study of self-gravitating systems from an equilibrium point of view. 

Then, a question naturally arises: Is what we learned with the above mentioned studies of regularized and confined 
models relevant to real self-gravitating systems, where no confining box is present? The aim of the present paper is 
to provide a first, preliminary answer to this question by studying the caloric curve of a basic, simple, yet effective 
model of globular clusters, the King model ^16] . This model describes a self-consistently confined system, without 
the need of any external box, but the stationary state is a non-thermal one. In particular, we shall consider the King 
model with a short-distance cutoff on the interactions and we shall discuss how such a cutoff affects the caloric curve. 

The paper is organized as follows. In Sec. [IT] we briefly review the known results on the equilbrium statistical 
mechanics of self-gravitating systems relevant to the understanding of our work. Section IIIII contains a discussion of 
the King model of star clusters as well as a calculation of its caloric curve. Section |IV] is devoted to the King model 
with a short-distance cutoff and contains the main results of the present work, especially concerning the caloric curve 
of that model and its dependence on the cutoff. Finally, Sec.|V]is devoted to some concluding remarks. 

II. EQUILIBRIUM STATISTICAL MECHANICS OF SELF-GRAVITATING SYSTEMS 

Let us consider a system of N classical point particles of mass m, mutually interacting via gravitational forces. The 
Hamiltonian, written as a function of the positions r^ and of the velocities v^ of the particles, is 

N N N 

■H(ri,...,rAr,vi,...,vjv) = y X! ^»^ ~ '^"^^ X! X! |r-~r I ' *'^'' 



where Vi = |vi| and G is Newton's gravitational constant. The equihbrium properties of such a system at constant 
energy E can be derived from the microcanonical (Bohzmann) entropy|43j 

S{E)=logoj{E) (2) 

where uj{E) is the density of states, defined by 

w(^)= /"<5[H(ri,...,rjv,vi,...,VAr)-^]dr, (3) 

with dT = rii^i '^^i '^^i- Equihbrium properties at fixed temperature T stem from the Helmhohz free energy 

F{T) = -T log Z{T) (4) 

where Z{T) is the canonical partition function given by 

Z{T) = j exp [-m (ri , . . . , r^, vi, . . . , v^)] dV , (5) 

where /3 = T^^ . As already discussed in Sec. IH the integrals in Eqs. ^ and ([S]) are well-defined only when the 
interaction potential is regularized at small length scales and the system is enclosed in a finite volume V . 

A. The isothermal sphere 

The analytical calculation of a;(ii^) or Z{T) for any system ([1]) with A^ > 2 is not feasible. The simplest approach to 
the equilibrium statistical mechanics of such a system when N is large is the mean-field approximation, which becomes 
exact in the A — >■ oo limit [3, la |lj| . The mean-field approach can even be pursued without the introduction of a 
short-distance regularization: in this case, the mean-field physical quantities remain finite, although their derivation 
from Eqs. ^ and @ is somehow ill-defined (see e.g. [8]). As we shall see in the following, however, neglecting the 
regularization does have a big effect on the physics. 

Let us now briefly recall an important example of the mean-field approach to the equilibrium of self-gravitating 
systems, the so-called "isothermal sphere", first studied in a seminal work by Antonov Tr]. This will also allow us to 
set up some notation which will be useful in the following. The system is enclosed in a spherical box of radius R and 
its state is described by a single-particle distribution function /(r,v) such that the total mass M is 



M = I f{r,v)drd\r, 



and its microcanonical entropy is a functional of / given by 



S[f] 



/(r,v) 
/o 



log 



/(r,v) 



/o 



drdw , 



(6) 



(7) 



where /q is a constant such that ///o is dimensionless, and otherwise arbitrary. Microcanonical equilibrium states 
correspond to maxima of the entropy. Here we have not introduced any short-range regularization of the potential, 
so that strictly speaking the entropy functional (O does not have a global maximum corresponding to a smooth 
distribution function /, because a singular, collapsed state would have a larger (actually infinite) entropy 0, [l4| . 
However, one can look for local maxima of the entropy and consider them as our equilibrium states, although they 
are only metastable states. Imposing the stationarity condition on the entropy one finds spherically symmetric 
distributions, i.e., /(r, v) ~ f{r,v), which depend only on the single-particle energy: 



f{r,v) = Cexp-^ -7 



(p{r) 



(8) 



where C is a normalization constant related to the total mass M, 7 is proportional to the inverse temperature and 
ip{r) is the mean-field gravitational potential which obeys the Poisson equation 



V^ip{r) = AnGgir) . 



(9) 



where g is the mass density 

g{r)^ I f{r,v)dv. (10) 



Finding the / corresponding to the local extrema of the entropy amounts then to solving the Poisson equation ([9|) , 
which becomes effectively one-dimensional due to the spherical symmetry, coupled to Eq. (j8]) via Eq. (fTO]). This 
cannot be done analytically but is an easy task on a computer. 

The collective behaviour in the microcanonical ensemble is conveniently encoded in the caloric curve, i.e., the 
relation between the temperature T and the energy E. To define these quantities we observe that the mean-field 
kinetic energy K is given by 

K = - v^f{r,v)drdv (11) 



and is related to the temperature T = (dS/dE) ^ via 



while the mean-field potential energy U is 



A-^f, (12) 



Cf = 5 / eWsJMdr, (13) 



so that E = K + U = ^ + U. Introducing dimensionless energy and temperature as 

" = CAP <'^' 

the caloric curve is obtained by plotting ^ as a function of e. The caloric curve of the isothermal sphere is reported 
in Fig. [T] The density contrast q{Q)/ q{R) grows along the curve, starting from q{Q)/ q{R) = 1 at e = -|-cx). At 
high energies ^ is approximately proportional to e and the systems is in a gas-like phase. By decreasing the energy, 
the temperature reaches a minimum value i?J„°,j — 0.397 and then begins to grow, so that the specific heat becomes 
negative. A minimum value of the energy e^^^^ ~ —0.335, corresponding to a density contrast g{0)/g{R) — 709, is 
eventually reached. The caloric curve then starts spiralling until a singular sphere is reached, with infinite density 
contrast. However, all the states with density contrast larger than 709 are unstable, because they are saddle points of 
the entropy rather than local maxima [^, [iJI ■ The part of the caloric curve corresponding to unstable states is plotted 
as a dashed line in Fig. [T] 

The presence of a minimal energy emin below which no equilibria are present is customarily interpreted as the 
"gravothermal catastrophe" : when e < ej^°j^ the system would collapse to a singular state. However, the presence of 
the gravothermal catastrophe is a consequence of the fact that we did not regularize the interactions at small length 
scales: as we shall see in the following, the introduction of a short-distance cutoff on the interaction potential removes 
the gravothermal catastrophe and stabilizes a low-energy phase. 

B. Confined models with a siiort-distance cutoff 

As already discussed in Sec. U a short-distance cutoff on the interactions must be introduced not only to make the 
density of states ([3]) and the partition function ([5]) well-defined [44], but also because it is physically motivated: no 
real system exists where the interaction potential is the classical gravitational one at all length scales. 

In the last decades, equilibrium statistical mechanics has been studied for many models of self-gravitating systems, 
defined in a confined volume and with small-scale regularization of the potential. Without entering the details, let us 
recall some of the latter models in order to underline their many common features, despite the models may be very 
different from each other. 

The isothermal sphere with cutoff, i.e., the model discussed above with the potential modified by the introduction 
of a small-scale cutoff, was first considered, in an approximate way, in 18]; more recently it has been fully studied 




FIG. 1: Caloric curve of the isothermal sphere. The red dashed line corresponds to unstable states. 



numerically, see e.g. [3, ISJ (also for 1/r" potentials with a 7^ 1) and [9|, ll3|. In these works the regularization of the 
potential is achieved via the so-called Plummer softening: the potential energy is given by 



with 



V(ri,...,rAr) 



Va{x) 



N N 
i—l j>i 



Gm^ 



v^ 



(16) 



(17) 



where a is the small-scale cutoff. The isothermal sphere with cutoff is a mean-field model, but the statistical mechanics 
of N particles mutually interacting with the potential energy given by Eq. (J16p has been studied also via microcanonical 
and canonical Monte Carlo simulations [1S|. A different kind of regularization has also been considered, i.e., a basis 
function regularization: the softening is achieved by truncating to n terms an expansion of the Newtonian potential 
in spherical Bessel functions [20| (a similar regularization was used in TV-particle dynamical simulations [21|). Also 
self-gravitating fermions have been studied, where the cutoff is not imposed but is a consequence of the Pauli exclusion 
principle J9l[l4|. 

As is customary in statistical mechanics, many simplified models have been proposed which try to grasp the essential 
physics by simplifying either the interactions or the geometry or both. A really minimalistic model was proposed by 
Thirring in the 1970's [23, where the only feature of the gravitational interactions that is retained is nonadditivity. 
Other popular models where the geometry is very simplified are the shell model [23| and the self-gravitating ring 
(SGR) 2^, 25i\ ; from the latter a minimalistic model of the same kind as Thirring's was recently derived ^2^] , showing 
that minimalistic models can be accurate also at a quantitative, and not only qualitative, level. Another simple model 
is the binary star model first introduced by Padmanabhan (see [8] and [14]), where only two particles are considered. 

All the above discussed models, although very different from each other, do share two main features: they are 
confined, in the sense that the configuration space is a compact manifold (either without boundary, as in the case of 
the SGR model, or with boundary when the systems is defined in R^ and confined in a box), and the potential is 
regularized at small length scales. The detailed properties of these models can be very different and also depend on 
the cutoff scale, and we refer the reader to the original works for these details. This notwithstanding, as anticipated 
in Sec. HI when the cutoff length is sufficiently small the caloric curve of all these models is qualitatively almost the 
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FIG. 2: Sketch of a typical caloric curve of a model of a regularized and confined self-gravitating system. 



same. At high energy there is a gas-hke phase, where the caloric curve is close to a straight Hne with positive slope. 
Lowering the energy the specific heat becomes negative. In some cases, the gas-like and the negative specific heat 
regions are separated by a microcanonical phase transition whose order does depend on the details of the model and 
on the specific short-distance regularization chosen. In other cases, a transition occurs within the negative specific 
heat region, i.e., there is a phase transition between two different negative specific heat phases [45J|. In these cases, the 
phase transition can be interpreted as the counterpart of the gravothermal catastrophe for models with short-range 
cutoff. The negative specific heat region, where temperature grows while decreasing the energy, is dominated by 
gravity; here the equilibrium state is clustered but typically a dense cluster is accompanied by a diffuse halo. At even 
lower energies the caloric curve reaches a maximum and the curve bends down, with a specific heat which becomes 
positive again and remains positive until the minimum energy is eventually reached. In this low-energy phase the 
physics is dominated by the regularization and the equilibrium state is strongly clustered. We note that, at variance 
with the isothermal sphere, the minimum energy of these models is the true minimum of their potential energy. A 
typical caloric curve of a model of a regularized and confined self-gravitating system is sketched in Fig. [21 

In the canonical ensemble, the intermediate negative specific heat phase is suppressed and there is a discontinuous 
phase transition between the gas-like and the low-energy collapsed phase. 

The general conclusion that can be drawn from these studies is that the introduction of a short-distance cutoff 
stabilizes a low-energy phase, which is the continuation of the negative specific heat region of the isothermal sphere. 
In this region the specific heat remains negative until eventually the regularization of the potential becomes dominant 
and temperature starts to decrease with energy until the minimum is reached: the smaller the cutoff length, the larger 
the negative specific heat region is. The gravothermal catastrophe in the microcanonical ensemble disappears. 

However, the presence of a high-energy gas-like phase, where the equilibrium state is nearly homogeneous, is a 
consequence of the confinement of the system. Hence we do not expect such phase to be observable in real self- 
gravitating system, except for some special cases where the system is indeed confined in a restricted volume by 
dynamical or other constraints I4g|. Clearly, also a sharp transition between the gas-like phase and the negative 
specific heat phase would not be observable, while a transition between different collapsed phases would not be ruled 
out. More generally it is unclear whether the "universal" collective behaviour observed in models of confined and 
regularized systems does indeed resemble the behaviour of real self-gravitating systems. In the next sections, and 
especially in Sec. |lVl we shall address this problem by studying the caloric curve of a model with small distance cutoff 
but no confining box. Such a model is a modification of a classic model of globular clusters, the dynamical King 
model [16,] : hence the next Section is devoted to a brief introduction to this model from the point of view of statistical 
mechanics. 



III. KING MODEL OF STAR CLUSTERS 

Globular clusters are systems of 10^ -^ 10^ stars, orbiting galaxies. Nearly 150 globular clusters orbiting the Milky 
Way have been discovered to date; around 500 globular clusters are estimated for the Andromeda galaxy M31 and 
roughly 10"' for the giant elliptical galaxy M87 in the Virgo cluster [gJl- Their shape is spherical to a very good 
degree of approximation and their size is finite: the limiting radius is referred to as "tidal radius" (denoted by rt) 
because for r > rt the tidal effects due to the host galaxy gravitational field dominate over the clusters's own gravity. 
Globular clusters are believed to be composed almost exclusively by stars, without relevant amounts of dust or gas 
(and dark matter too). This makes globular clusters almost ideal self-gravitating systems. Moreover, they are very 
old: the age of globular clusters is typically larger than 10"'^° years [27|, which is also the order of magnitude of the 
coUisional relaxation time as estimated via the classical Chandrasekhar formula, so that it is reasonable that they 
have undergone some kind of relaxation. 

Globular clusters are therefore the natural choice as examples of real self-gravitating systems that could be ap- 
proached by statistical mechanics. However, standard thermal equilibrium can not be present in these systems, since 
they are not enclosed in a box: their velocity distribution will not be a pure Maxwellian one. However, it is reasonable 
to assume the velocity distribution is as close to a Maxwellian as possible. Indeed, the idea at the basis of the King 
model is a modification of the isothermal sphere discussed in Sec. lH AK see e.g. [l[). The state of the system is defined, 
as in the isothermal sphere, by a spherically symmetric single-particle distribution function /(r, v) such that the total 
mass M is given by Eq. ([S]). The velocity distribution is obtained by "lowering a Maxwellian^V, such that no particle 
can have a velocity larger than the escape velocity Ve{r). The King distribution function is [l6| 



J(r,v) = < I J 

lo iiv^>vl{r), 



(18) 



where C is proportional to the total mass M and 7 is an inverse potential scale which is no longer directly proportional 
to the inverse temperature as the 7 of the isothermal sphere was, since the velocity distribution function is no longer 
Gaussian. Let us choose energies such that ip{rt) — 0, with the tidal radius rt defined as the radius such that the 
density g, still given by Eq. (jlOp . vanishes; then Wg(r) = —2(p{r). Given these definitions, everything goes on as 
in the case of the isothermal sphere; since we have not introduced any short-distance cutoff on the iteractions (we 
shall consider this case in Sec. lIVp the potential Lp{r) obeys the Poisson equation ^ and we get an effectively one- 
dimensional differential equation whose solution gives the explicit /. The difference with the case of the isothermal 
sphere is that the / thus obtained is no longer a stationary point of the microcanonical entropy [43|. In practice one 
chooses a value for the central potential energy (pQ or for the central density go = p(0) and integrates the differential 
equation to get the density profile g{r); we refer the reader to the original paper [1^ and to [ll, |23| for the details. If 
rt — >■ cx), then the King model becomes identical to an isothermal sphere with i? — > cxd and M —> 00. 

Assuming that "light follows mass", i.e., that the luminosity density is proportional to the mass density, one can 
project g{r) onto the plane of the sky [l|, |lg to compare the result with photometric observations. The resulting 
distributions provide very good fits [48] to the observed densities of nearly 80% of the Milky way globular clusters 
[la [29I . l3g | . Also the internal kinematics (the velocity dispersion profile) of the globular clusters is reasonably well 
described by the King model [29|, [3l|, even if there may be significant deviations (see e.g. [S^). Although the recent 
availability of higher- quality data has shown that the King model does not provide an optimal description of globular 
clusters, it is still considered the best physically motivated model to date (see e.g. [31|) and it remains the standard 
data interpretation and classification tool in this field |28| . 

A. Caloric curve 

The kinetic and potential energy of the King model can be calculated as in the case of the isothermal sphere. Kinetic 
energy is still given by Eq. (jlip . while potential energy is given by the same expression as in Eq. (|13p . provided the 
potential is not the potential f{r) entering Eq. (jl8|) but the potential vanishing at infinity, ipocif) = fif) — ^(oo) = 
ip{r) — GM/rt; the potential energy of the King model can thus be written as 

t/ = h27r / g(r) ip(r) r^ dr . (19) 

'2rt Jo 

Then, the temperature T may still be defined as in Eq. (J12p . although we do no longer have a microcanonical entropy 
such that T — (dS/dE)^^ . Hence a "microcanonical caloric curve" is still given by plotting temperature versus 
energy, although strictly speaking there is no equilibrium microcanonical ensemble; the constraint of constant energy. 



however, justifies the analogy. It remains to define the dimensionless energy e and temperature i?, and in close analogy 
to the case of the isothermal sphere we set 

- '^'^ (20) 



GAP 



where the tidal radius rt plays the role of the radius R of the box of the isothermal sphere. Using r^ as unit of lengths 
is not customary in the context of the King model, where a length scale derived from dimensional analysis of the 
Poisson equation is typically used [16ll3l|: however, our choice simplifies the comparison between the King model and 
the isothermal sphere, because one is somehow using the same unit of length in the two cases. 

At variance with the isothermal sphere, the shape of the caloric curve 'd{e) of the King model can be predicted 
without performing any calculation: since the King model gives a stationary distribution function and involves only 
pure gravitational interactions, without external containers or short-distance cutoffs, the virial theorem explicitly 
gives the relation between temperature and energy as 

^ = -^e • (22) 

However, the virial theorem alone does not tell anything on the boundaries (if any) of the domain of the caloric curve. 
From Eq. (J19p we have that the dimensionless potential energy of the King model 



GA'P 



(23) 



is surely bounded above by u = ^1/2, but the calculations show that the real upper bound is smaller, and turns out 
to be Wmax — —1-20. Hence, using again the virial relation K = —U/2 we have that e < e^™^ ~ —0.60 in the King 
model. Performing the calculations shows that also a lower bound on energy exists, i.e., e G [£min®,emax'] where 



^min — ^-ro, 



^^!Sf ^ -0-60. 



(24) 



The latter result means that no King distribution (|18p exists such that the dimensionless mean-field energy is larger 
than ejj"ja5f or smaller than s^^'^. Due to the virial relation (|22]) also the temperature is bounded, i.e., ■& E [^mln^' i^^'"^^ 



with 



^Kmg ^ _ 2 King ^0 40 

^miii 3^max — "-/.-^v^ , 

,QKing ^ _ 2 King ^ -, ,„ 

''max a'-min — j-.t:^ . 



(25) 



This is at variance with the isothermal sphere, where the temperature is unbounded above due to the existence 
of the gas phase; but if we consider only the negative specific heat region of the isothermal sphere, also in that 
case the dimensionless temperature is bounded from above an below; incidentally, the two lower bounds are very 
similar. Hence, the qualitative appearance of the caloric curve of the King model is not very different from that of the 
isothermal sphere, provided one removes in the latter the gas-like phase due to the presence of the container and one 
takes into account the virial theorem constraint; the King model minimal energy e^^^ ~ —2.13, is, however, sizably 
smaller than the corresponding value of the isothermal sphere e[^°jj ~ —0.335. 

The caloric curve of the King model is reported in Fig. |3l the straight line is the virial relation (f22|) . while 
the symbols are (e, 1}) values calculated using King distribution functions. The accumulation of points close to 
(£,i9) = {e*,-d*) ~ (—1.65, 1.1) in Fig. [3] is a remarkable effect. Indeed, by calculating King distribution functions for 
increasing central densities qq — g{r — 0), the corresponding values of e decrease until e^-^* is reached. Continuing 
to increase go results in an oscillating pattern of e values centered around e = £*, as in a "collapsed spiral". The 
value e = e* is asymptotically reached for qq -^ 00. This behaviour closely reminds of the spiral behaviour of the 
caloric curve of the isothermal sphere: here the spiral collapses onto a straight line due to the constraint (|22|) . To 
show that it is really a collapsed spiral, we plot not "d but a related quantity as a function of e. First, we define a 
"central temperature" Tq by 

Ko='-^ (26) 



=& 




FIG. 3: Caloric curve of the King model. The solid line is the i9 = 
calculated using King distribution functions. 



je virial law and blue circles are the values of ■& and e 



where Kq is the kinetic energy the system would have if v^ were constant and equal to that in the center [4| 



v\r = Q), 



Kn = -Mv, 



Then, in Fig. |4] we plot the dimensionless central temperature t?o given by 



i?o = 



GNP 



(27) 



(28) 



as a function of e. Having removed the constraint that the values must lie on a straight line, the spiral pattern opens 
up. 

Although not expressed in terms of a caloric curve, the presence of a minimal value for the energy was already 
noticed in King's original paper |16| and interpreted as a stability limit for the model, suggesting that the model 
should be modified to accommodate more collapsed density profiles (the issue of the stability limits of King and other 
similar models was further elaborated in [33l| and [3J|). From the above discussion it is apparent that the presence of 
a minimal energy e^^ is nothing but the analogue in the King model context of the gravothermal catastrophe of the 
isothermal sphere. Hence, in view of what recalled in Sec. IHBl we suggest that the addition of a small-scale cutoff to 
the interactions is likely to be a natural and physically motivated modification of the King model which stabilizes a 
low-energy phase where more collapsed density profiles are allowed. 

The next Section is devoted to discussing the King model with a short-distance cutoff. 



IV. KING MODEL WITH A SHORT-DISTANCE CUTOFF 



Let us now define the King model with short-distance cutoff on the interactions. We assume that the distribution 
function is the King one (|18|) . We choose a Plummer short-distance regularization of the interactions, given by Eq. 
(|17p . The regularized gravitational potential i/'oo(r) generated by a mass density g can thus be written as 



V'oo(r) = -G / 



:dr': 



(29) 



r -|- a^ 
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FIG. 4: Central dimensionless temperature ^o of the King model as a function of the dimensionless energy e (blue solid line), 
compared to the caloric curve i9(e) (black dash-dotted line) already reported in Fig. (3] 



as before, we denote by (poo the gravitational potential that vanishes at r — > oo, to distinguish it from the potential 
Lp entering the expression (|18p of the distribution function, that vanishes at rt- Since our density is spherically 
symmetric, also the potential must depend only on r. Indeed, for p(r) = Q{r) the angular integrations in Eq. ([29)1 can 
be performed to yield 



Vooir) = - 



2ttG 



r' g{r') y {r' + r)-^ + a^ — \/{r' — r)^ + a? dr' , 



(30) 



where we have taken into account that g{r) = if r > r^. The potential (p{r) entering Eq. (1181) can be obtained from 
Eq. ((30)) by subtracting ipooift) and reads as 



(p{r) = ~2ttG / r' g{r') 
Jo 



■y/(r' + r)^ + a^ — ■\/(r' — r)^ + a^ ■\/(r' + rt)^ + a^ — ■\/(r' — rt)^ + i 



n 



dr' . (31) 



Inspection of Eq. pip shows that (p{rt) = 0. 

At variance with the standard King model, our potential does not obey a Poisson or Poisson-like equation. Hence 
to find the density profile g{r) explicitly we have to solve self-consistently the equation giving the density as a function 
of the potential, that is Eq. (fTO)) with / given by Eq. (fTSt , and the equation (|3Tt defining the potential as a function 
of the density. To this end it is convenient to define a dimensionless potential W as 



W{r) = -2-fip{r) 



(32) 



and an analogous Woo defined as above with ipoo in place of (p. In Eq. (15^ 7 is the inverse potential scale entering 
the definition (fT5)) of the distribution function. By definition, W > ii r < rt with W{rt) = 0. We then integrate Eq. 
(ITU)) with / given by Eq. ((TH]) to obtain 



g{r) 



47rC 



$ [W{r)] 



where 



$(iy) = 2e 



w 



w 



'V „i 



rj drj 



(33) 



(34) 
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The normalization constant C can be expressed as a function of the total mass M, 

~ 167r2 /J' $ [W{r)] r^dr ' ^ ' 

Since W{rt) = 0, Eq. p3|) correctly implies e(rj) = 0. The two equations to be solved self-consistently are thus Eq. 
([55)1 and Eq. ([5^ . with ip{r) given by Eq. (pij) . It is convenient to rewrite them using dimensionless quantities. Let 
us put rt as unit of length, M as unit of mass and GM"^ /rt as unit of energy, in order to use the same units as in the 
previous Section. We then define a dimensionless radius x as 



a dimensionless mass density V' 



a dimensionless inverse potential scale p 



X 


r 

n 


i;-- 


grl 
M ' 


■ 2 

7 = 


7GM 



(36) 



(37) 



(38) 
n 

and a dimensionless cutoff length a 

a = - . (39) 

n 

Using these quantities the system of dimensionless equations to be numerically solved self-consistently can be written 
as 

^(.) ^ -^^^ (40) 

W{x) = W^{x)-WU^) (41) 

WU^) ^ 1^ /■ y^(y) \^{^y + n,y + a^^^{y-xY + aA dy (42) 

•^ Jo L -I 

with $ given by Eq. (I34p . The above equations have to be solved iteratively, starting from an initial condition for the 
dimensionless density ipo{x)- During the iterative procedure the constraints that the density and potential vanish at 
a; = 1, i.e., V(l) = W{1) = 0, and 

47r / V(a;) x^dx = 1 , (43) 

Jo 



i.e., that the total dimensionless mass is equal to one, must be satisfied. Due to the form of Eqs. (|40)) . pT|) and p2)) 
these constraints are automatically satisfied during the iterative procedure once an initial condition ipo (x) compatible 
with the constraints is chosen. 

In practice, after fixing a value for a, we choose a value Wq and set j^ such that W{0) — Wq; then we set ip to an 
initial form[50| V'o(2;) which satisfies the constraints. Then we compute W from Eqs. (|4T|) and p2|) and we insert the 
function W{x) thus obtained into Eq. (j40p . finding a new dimensionless density ipix) which is inserted again into Eq. 
(H^ . We iterate the procedure until the integrated difference between densities at two subsequent steps is smaller 
than a prescribed threshold i51j. This yields the desired dimensionless density profile corresponding to the chosen 
value of Wo- Then we choose another value for the control parameter Wq and repeat the procedure, starting from the 
ip{x) found at the previous value of Wq: we also checked that the result does not change if we start from ^'0(2:)- 

A. Caloric curve 

Given a density profile ip{x) and the corresponding potential profile W{x) we can calculate the dimensionless kinetic 
and potential energy, from which we obtain the dimensionless temperature 1? and energy e. 
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The dimensional potential energy U is 

l-rt 

U^2n r'^ if^ir) g{r) dr , (44) 

Jo 

so that the dimensionless potential energy u = Urt/{GA''P) is given by 

u = -— WUx)i^{x)x'dx. (45) 

J^ Jo 

The dimensional kinetic energy K can be written as 

K = 8tt^ drr^ dvv^f{r,v); (46) 

Jo JO 

with the change of variable 77 = 71;^ and one partial integration we get 

K = ,,^ ■^'\ ^^ ^ / dr r" e^f"^) / ^''^ e"" dry , (47) 



hence the dimensionless kinetic energy k = Krt/{GM ) is 



o i-l i-W{x) 

' dxx^e"^^^^ / 77^/2 e-^dr?. (48) 



mf j^<^[W{x)]x'^dxJo Jo 

The dimensionless temperature d = Trt/{GM'^) is then given by 

^=1^, (49) 

and the total dimensionless energy s is 

e = K + u . (50) 

The caloric curve i? vs. e of the King model with short-distance cutoff obtained as explained above is reported in 
Fig. [5] for a dimensionless squared cutoff length a^ = 10^'^ and in Fig. [6] for a smaller value of the cutoff, a^ = 10~^. 
Two facts are apparent: (i) the small-scale cutoff stabilizes a low energy phase extending until energies well below 
the minimal energy e„-"® ~ —2.13, i.e., until the minimum value 

o.y mill ' " 

^" — ,," = (r^-\) 

^min '^min 9 ' \^ J 

corresponding to the potential energy of a system where all the particles are in the same point, so that such a low- 
energy energy phase gets larger as the cutoff decreases; (ii) the shape of the caloric curve is very similar to that of 
confined and regularized models, discussed in Sec. IIIBI and sketched in Fig. [2l but for the absence of the gas-like 
behaviour at high energies, which was allowed in those models by the presence of a confining box. 

In Fig. [7] a close-up of the caloric curves with cutoff (for the two cutoff lengths a^ — 10"'^ and a^ = 10^^ of Figs. [S] 
and [HI respectively) and without cutoff is shown in energy interval containing the minimal energy of the King model 
without cutoff, £jj,in^ — —2.13. We observe that already for a^ = lO^'"' the caloric curve with cutoff is practically the 
same as that without cutoff, where the latter exists, but extends well below the latter as already shown before. 

These results give a first and positive answer to the question we put at the beginning of the present work, i.e., whether 
the behaviour of regularized and confined models of self-gravitating systems might be relevant to real systems where 
neither confinement nor proper thermal equilibrium states exist. Indeed, by studying a short-distance regularized 
version of an observationally probed model, the King model, we have shown that the caloric curve of such a model is 
qualitatively the same as that of regularized and confined model, the only qualitative difference being the absence, in 
the King model with cutoff, of a gas-like phase at high energies. This difference was expected from the outset since 
the gas-like phase is allowed only due to the confinement. Moreover, we have shown that already for moderate values 
of the cutoff (a^ = 10^^, see Fig.[71) the caloric curve of the King model with cutoff is virtually indistinguishable from 
that of the original King model without cutoff, for values of energies where both curves exist. 

In the next two subsections we study the shape of the density profiles obtained with the King model with cutoff 
and we address the problem of understanding what happens when the cutoff gets smaller than the values considered 
up to now. 
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FIG. 5: Caloric curve of the King model with a short-distance cutoff. The dimensionless squared cutoff length is a^ — 10 



B. Density profiles 



We can recognize three different regions in the caloric curve of the King model with cutoff. First, a high-energy 
region, allowed also for the model without cutoff, for values of e > ^min^- Second, an intermediate-energy region, which 
is the continuation of the caloric curve without cutoff for e < ^min^ ^^^ where the specific heat remains negative. 
Third, a low-energy region where the specific heat becomes positive and is the cutoff-dominated region. 

It is interesting to have a look at the density profiles 'ip{x) corresponding to the three different regions. Examples 
of density profiles are reported in Figs. [8] and [9] for a^ — \Q^^ and a^ — 10^^, respectively. Note that the density 
profiles are plotted using log-log scales. The main feature emerging from these plots is that there is a qualitative 
change in the shape of the density profile passing from the first to the second and third regions of energy. In the 
high-energy region, where for a^ — 10^^ the caloric curves with cutoff and without cutoff are nearly indistinguishable, 
the density profile is, as expected, very close to the density profile of the King model without cutoff 52], i.e., it starts 
with a fiattened core and then decays without infiections. When we lower the energy an infiection point in the profile 
appears and a central core develops, whose density grows as the energy is lowered. For the smaller of the two cutoffs, 
the change between the two kinds of density profiles (without and with the central core) is much more abrupt. 

The interest of this phenomenon is twofold. First, it suggests that the physics in the energy region that was 
not allowed in the model without cutoff is not a simple continuation of the previous one and that at smaller cutoff 
the change between the two kinds of density profiles might correspond to a phase-transition-like phenomenon; this 
possibility will be addressed in the next Section. Second, it has also a phenomenological interest. We have mentioned 
that the King model provides good fits to nearly 80% of the Milky Way globular clusters. Most of the remaining 20% 
are referred to as "(post-)core-collapsed" clusters 35]. It is believed that, after having undergone a gravothermal-like 
collapse in the central region which has been stopped by dynamical effects like the formation of binary stars, these 
clusters have settled in a near-equilibrium state which cannot be described by the King model because their density 
profile is too much concentrated [36] . Many of these clusters have surface brightness profiles with a central core, 
qualitatively similar |53| to those we obtain with our model with cutoff. A detailed comparison with observations is 
outside the scope of the present work, and we are not claiming our model with cutoff is able to provide a good fit for the 
density profiles of core-collapsed clusters; however, the qualitative similarity between observed and theoretical profiles 
suggests that the introduction of a short-distance cutoff in the interactions might be a possible way to construct 
simple models able to describe also these clusters [54]. 
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FIG. 6: As in Fig. [S] with a^ = 10" 




FIG. 7: Caloric curves of the King model with a short-distance cutofT (blue dashed line: a^ = 10 ^, red dot-dashed line: 
o? = 10~^) compared with the caloric curve of the King model without cutoff (black solid line). The dotted line is the 
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FIG. 8: Dimensionless density profile tp{x) as a function of the dimensionless radius x for the King model with short-distance 
cutoff, with a^ = W~^. Black solid curve: e = —1.27, belonging to the energy range in common with the King model without 
cutoff; blue dashed curve: e = —3.91, belonging to the intermediate energy region with negative specific heat; red dot-dashed 
curve: e — —9.02, belonging to the lowest energy region with positive specific heat. The dotted vertical line marks the cutoff 
length scale a. 
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FIG. 9: As in Fig. [S] with a^ = 10 ^ and energies: e = —1.37 (black solid curve); e — —7.19 (blue dashed curve); e 
(red dot-dashed curve). 
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FIG. 10: Comparison between the caloric curve of the King model with a short-distance cutoflf with o? — 7.5 x 10 '^ (blue 
dot-dashed line) and the caloric curve of the King model without cutoff (black solid line) . The dotted line is the continuation 
of the virial law i9 = — |e for energies below e^^^ — —2.13. A change in the slope of the caloric curve with cutoff is apparent 
for energies between 2 and 3. 



Phase transition for small values of the cutoff? 



The analysis reported in Sec. lIVBJ has shown that in the low energy regime, i.e., when e < e^^^^, the density profiles 
of the King model with short-distance cutoff are qualitatively different from those at higher energies, the latter being 
essentially equal to those of the King model without cutoff. This suggests that for small a's there might be a sharp 



crossover or even a singularity in the caloric curve for e ~ e 



King 



This would be analogous to the phase transition 



between different states with negative specific heat that is observed e.g. in the isothermal sphere with short-distance 
cutoff |14| , for small values of the cutoff. 

There is another reason to expect that something special happens to the caloric curve close to £ ~ £min^' f*^^ small 
values of a. The caloric curve of the King model with cutoff must become equal to that of the King model without 
cutoff when a — >■ 0. We have already seen that, for energies in the domain of the caloric curve without cutoff, the 
caloric curve with cutoff becomes practically equal to that without cutoff already for a^ = 10~^, so that we expect 
this to be true for any a^ < 10"^. However, for finite a the domain of the caloric curve with cutoff appears much 
larger than that without cutoff, and getting larger as a gets smaller; but when a — > also the domains of the caloric 
curves must become the same, so that all the low energy region e < s^^^ ought to disappear. There could be, in 
principle, (at least) two different scenarios leading to that. In the first scenario the limit a — > is a singular limit: the 
whole low-energy region exists for any a > and abruptly disappears when a — >■ 0. In this case one could also expect 
that for sufhciently small a a phenomenon analogous to an equilibrium phase transition exists, where the caloric 
curve develops a singularity separating the low-energy phase from the high-energy phase whose domain contains (and 
when a — 7> coincides with) that of the King model without cutoff. In the second scenario, for sufficiently small a a 
forbidden region in energy eiow(o) < £ < £high(Q^) opens up, where no solution of the self-consistent equations (j40l) . 
(|4T]) and ((42]) exists. The extrema of the forbidden region in energy would be such that lima^o £iow(a) = Emin ^^"^ 
lim^^o £high(«) — ^min^' ^^ ^^^* although the minimum energy decreases with a, also the effective width of the low 
energy phase decreases and eventually shrinks to a point as a — ?► 0. 

The results we have obtained so far do not allow to discriminate between the two scenarios yet, although the phase- 
transition-like scenario seems favoured. The first result suggesting a phase transition might occur is shown in Fig. [TUl 
where we plot the caloric curve for a^ = 7.5 x 10~^, in a slightly wider energy window than in Fig. [71 and a change 
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FIG. 11: Dimensionless energy £ as a function of the control parameter Wo (dimensionless central potential) for the King model 



with a short-distance cutoff. Blue dot-dashed curve: 



10 ; red solid curve: a 



converged (non-oscillating) solutions; red dashed and dotted curves: a = 5 x 10 
lower bounds of oscillating solutions. 



= 5 X 10 , values of e corresponding to 
values of e corresponding to upper and 



of slope in 'd{e) at energies between 2 and 3 is apparent, which reminds of a caloric curve near a continuous phase 
transition. This suggests that for slightly smaller a the caloric curve may develop a kink and, for even smaller a, a 
discontinuity like in a first-order phase transition. However, upon slightly decreasing the cutoff length to a^ = 5x 10~®, 
we were not able to find a reliable convergence to a self-consistent solution for values of the control parameter such 
that 9.67 < Wq < 10.38. When trying to solve the self-consistent equations for these values of Wq, the solution did 
not converge yet also after a number of iterations orders of magnitude larger than the number of iterations leading to 
convergence in the rest of the curve. In these cases the kinetic and potential energies oscillate between two different 
values. In Fig. [TT] the values of e are plotted as a function of the control parameter Wq, for two choices of the cutoff, 
a^ = 10~^ and a^ = 5 x 10^^. In the latter case, in the region where no reliable convergence was found, upper and 
lower bounds of the oscillations are reported. This behaviour is compatible with both scenarios, because oscillating 
solutions may indicate either true absence of self-consistent solutions or a hysteresis-like phenomenon, due to the 
existence of metastable states, which in turn might suggest a first-order-transition-like behaviour. Indeed, even if no 
self-consistent solutions exists for Wq € [W'o°^, Wg '* ], this does not automatically mean that a forbidden region in e 
exists; this would require also that e^Wq""^) ^ ^(Wq '* )' while the shape of the curve reported in Fig. [TT] suggests that 
e(W^o°'") = £(M/(5"^'') is not ruled out. The caloric curve ^(e) for a^ = 5 x 10~^ is plotted in Fig.fTU sohd fines refer to 
converged solutions while dotted lines are upper and lower bounds of oscillating solutions. Although a small window 
in energy where no converged solutions exist is present (—3.120 < e < —2.715), bounds of oscillating solutions lie 
on continuations of the converged solutions and resemble metastable states typically observed close to discontinuous 
phase transitions. Indeed, finding the converged solution for energies where also oscillating solutions exists required a 
fine sampling of Wq values as well as a very large number of iterations: with the typical number of iterations leading 
to convergence in other regions of the curve, one would see oscillations. Therefore, at least for a^ — 5 x 10~^, the 
presence of the small energy region without any converged solutions might well be a purely numerical problem due 
to the very slow convergence of our numerical method in this region, which could perhaps be solved employing a 
much larger number of iterations [55|. Hence the phase-transition-like scenario seems favoured. Going to even smaller 
values of a yields a very similar behaviour, but for that the Wq range where solutions oscillate gets larger. 

The overall shape of the caloric curve for a^ = 5 x 10~^ is very similar to that shown in Fig. [5] so that we do not 
show it. Some density profiles obtained with a^ = 5 x 10^^ are plotted in Fig. [T31 We observe that the qualitative 
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FIG. 12: Caloric curve of the King model with a short-distance cutoff, in an anergy range similar to that of Fig. 1101 The 
dimensionless squared cutoff length is q^ = 5 x 10~®. The red and blue solid lines are converged values, while the red and blue 
dotted lines refer to the lower and upper bounds of oscillating solutions. 



features of the profiles are the same as for larger cutoffs, 
change the shape of the density profiles. 



Hence going to smaller values of a does not qualitatively 



CONCLUDING REMARKS 



The main motivation of the present work was to try to understand whether the behaviour of regularized and confined 
models of self-gravitating systems previously studied in equilibrium statistical mechanics could be considered relevant 
to real self-gravitating systems like globular star clusters, where there is no confinement so that the stationary velocity 
distribution is not the thermal one. To this end, we have studied a modified King model of a globular cluster, with 
a short-distance cutoff on the interactions. We have shown that such a model does behave like the regularized and 
confined ones, especially for what concerns the caloric curve, but for the absence of a high-energy gas-like phase which 
can only be induced by confinement. The caloric curve and density profiles of our model become equal to those of 
the observationally probed King model without cutoff for cutoff lengths such that a^ < 10^^, for energy values larger 
than the minimum energy allowed for the King model without cutoff. 

Our results thus provide a first, positive answer to the question that motivated our work. 

There are still many open points. First of all, the issue of the presence and of the nature of a transitional phenomenon 
for small cutoffs, discussed in Sec. lIVCl should be clarified. To this end it might probably be useful to consider different 
regularization schemes allowing for numerically easier treatments (possibly in a differential form). Moreover, although 
the model studied in the present work does not have a proper thermal equilibrium, one might construct an effective 
entropy in close analogy to Eq. ([7]) and try to use it to study the stability of the solutions. 

Another open problem concerns the physical interpretation of the small-scale cutoff. In the present work, as in 
previous works recalled in Sec. Ill Bl the short-distance cutoff has been taken as an external parameter, whose presence 
is required by physical considerations and whose numerical values should be sufficiently small to yield results not too 
far from those of the King model without cutoff at large energies and are otherwise arbitrary. However, physical 
considerations should also provide reasonable bounds for the cutoff scale or even estimates of its numerical value. 
Clearly the cutoff scale must be bounded below by the size of the individual stars; an order of magnitude for an 
upper bound could be given by the average interstellar distance. Translated into our dimensionless units this means 
-[^Q-9 < a < 5 X 10~^, so that all the cutoffs we have considered above fit into this range. Theoretical considerations 
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5 X 10 '^ and energies: e = —1.37 (black solid curve); e = —17.65 (blue dashed curve) 



suggest that the cutoff scale is more likely to be close to the upper bound than to the lower bound: since we describe 
the system with a single-particle distribution function, the cutoff scale should be of the same order of magnitude as 
the interparticle distance in order to be consistent with neglecting correlation effects (see e.g. [33|)- On the other 
hand, the cutoff might be considered as an effective implementation of dynamical effects that stabilize an otherwise 
unstable highly-concentrated distribution of mass, like the formation of binary stars: reasonable values of the cutoff 
may thus be suggested by analyzing such effects. 

As a final remark we recall that the King model with cutoff allows for energies much smaller than those of the King 
model without cutoff, corresponding to density profiles with a central core. From a phenomenological point of view, 
this suggests that also core-collapsed globular clusters not well described by the King model might be accommodated 
in a description with cutoff. A detailed comparison between the model predictions and observed density profiles is 
left for future work: in case the density profiles with cutoff turn out to yield reasonable descriptions of core-collapsed 
clusters, these clusters might be interpreted as the manifestation of a phase transition between different structures 
described by the same model, instead of as something requiring a completely different description. 
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